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ABS fRA aT Comparative modeling meth- 
ods arc described that can be used to construct 
a tnrec-dimensional model structure of a new 
protein trom l^nowledge of its sequence and of 
the experimemtal structures and sequences of 
other mertibei-s of its homology famUy. The 
methods are illustrated with the mammalian 
serine protease family, for which seven experi- 
mental structures have been reported in the lit- 
erature, and the sequences for over 35 different 
protein members of the family are available. 
The strategy for modeling these proteins is pre- 
sented, and criteria are developed for determin- 
ing and assigning the reliability of the modeled 
structure. Criteria are described that are spe- 
cially designed to help detect cases in which it 
is likely that the local structure diverges signif- 
icantly from the usual conformation of the fam- 
ily. 

Key words: protein structure, computer model- 
ing, structure prediction, sequence 
homology, structure homology 

INTRODUCTION 

It has been apparent for close to two decades that 
proteins from very different sources and sometimes 
with rather diverse functions can have homologous 
sequences and consequently very similar three- 
dimensional structures.^ This fact has been the ba- 
sis for the development of comparative modeling 
methods,'^"^ which permit extrapolation from the ex- 
perimentally determined structure for one or more 
members of a homologous family to a new member of 
this family whose sequence has been determined but 
whose structure is as yet unknown. A number of 
factors combine to increase greatly the application 
of comparative modeling techniques today. The 
large number of protein stnictures^ and the explod- 
ing number of protein sequences^ that are being re- 
ported in the literature and that are readUy ava.^ 
able in computerized databases provide the bas^^ 
structural and sequence data needed to apply the 
mlthld The proteins in which we are interested are 
method. Tbe PJ^ small quantities, too little for 

often available in o ^^^^^ 

structural ^^"^J^^^^^ and expressed. Even when this 
S^r e^rtls deemed worthwhile and is initiated, 
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the comparative modeling studies can be performed 
in the meantime, providing an approximate view of 
the structure of the molecule until sufficient protein 
can be, obtained and the experimental structure can 
be determined. Such a model structure can be very 
useful in the interim to help plan and interpret 
biochemical and mutagenesis-*-^ experiments, to 
probe functional properties, and improve our un- 
derstanding of ligand or substrate binding. ^^"^^ 

The use of comparative methods involves extrap- 
olation from one or more known structures to con- 
struct a new structure. It is important to consider 
how accurately this function can be performed in 
assessing whether it is a useful and worthwhile ex- 
ercise. Certainly, comparisons that have been pub- 
lished between predicted model structures and the 
experimental struct xires determined later ^^"^^ indi- 
cate that the modeled structure is not likely to be 
completely accurate. This paper reports the methods 
that we have developed to perform comparative 
modeling and the criteria that we use to estimate 
the reliability of various parts of the structure and 
especially to identify potential problem areas that 
are likely to be particularly difficult to predict cor- 
rectly. 

To describe properly and illustrate the modeling 
techniques, we will use the mammalian serine pro- 
tease family. Members of this family are ubiqui- 
tous in nature; they are present all along the evolu- 
tionary pathway from bacteria to humans. They 
play important roles in a wide variety of body func- 
tions, including blood coagulation, fibrinolysis, com- 
plement activation, fertilization, and digestion. 
Therefore, modeling the various members of this 
family can provide valuable new information about 
divers critical biological functions of the body. 

Comparative modeHng works best when there are 
several experimental structures and known se- 
quences. Experimental structures for seven differ- 
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TABLE I. Experimentally Known Three- 
Dimenslonal Structures of Serine Proteases 







Resolution 


Refer- 


Protein 


Source 


(A) 


ence 


Chymotrypsin 


Bovine 


1.8 


24* 


Trypsin 


Bovine 


1.8 


25* 


Elastase 


Porcine 


1.8 


26 


Kallikrein 


Porcine 


1.8 


27 


Mast cell protease 


Rat 


1.8 


28 


S. griseus trypsin 


S. griseus 


1.7 


23 


Tonin 


Rat 


1.8 


29 



*Structures have been determined by several groups for each of 
these proteins. The reference shown corresponds to the coordi- 
nates used in this work. 



ent serine proteases* that can be considered mem- 
bers of this mammalian family have been reported 
in the literature and deposited in the Brookhaven 
Protein Structure Database.® These include chy- 
motrypsin,^'* trypsin,^^ elastase,^® kallikrein,^'' rat 
mast cell,^® Streptomyces griseus trjrpsin-like 
protein,^^ and tonin^® (see Table I). Several of these 
have been solved independently in more than one 
laboratory either in the same form or in alternate 
forms of the molecule. ^^"^^ In addition, structures 
have been reported for the precursor, zymogen, form 
of chymotrypsin^^*^"* and trypsin,^^-^® To comple- 
ment these structures, amino acid sequences for 
over 35 different serine proteases have been re- 
ported (Table II), not counting species variations. 
Thus this family presents an ideal system for devel- 
oping comparative modeling methods and at the 
same time applying them to proteins involved in im- 
portant and interesting biological functions. 

COMPARATIVE MODELING METHOD: THE 
"SPARE PARTS" ALGORITHM 

Comparative modeling requires extrapolation 
from known structures to produce a model of an un- 
known structure. Our ability to extrapolate accu- 
rately is, unfortunately, greatly limited by our still 
rudimentary knowledge and understanding of pro- 
tein structure and energetics. Consequently, the 
techniques that are used to extrapolate influence 
the resultant model critically and may introduce 
considerable error. The methods that we have devel- 
oped attempt to systematize the modeling process by 
combining all the known structures and sequences 
to help improve the accuracy of the extrapolation. 



:;^v ?*Structures have also been reported for several serine pro- 

;^^^K:>VA:^fmit^^s^s froni bacterial sources,^'^ Although these are cer- 



Although these are cer- 



%0St^^'i^Sh\y members of this family, they differ sufficiently in.Uieir 
V^lJ^tS 'SI's^^^ that comparative modeling methpik. are not^de^^^ 



Analysis of the Structural Properties 
of the Family 

The first steps in the modeling process are illus- 
trated schematically in Figure 1. Let us assume that 
we have experimentally determined three-dimen- 
sional structures for several members of the homol- 
ogous family of interest, e.g., structures "A," "B," 
and "C" in Figure la. The known structures are su- 
perimposed in three dimensions to obtain a maximal 
overlap of the structures (Fig. lb). Performing this 
superposition of the structures is sometimes not a 
trivial procedure. This is true because we want to 
overlap the molecules based on the parts that are 
the same and ignore the parts that are different. The 
more different the molecules are from each other, 
the harder it is to find a unique best overlap of the 
common features of the structures. Several pro- 
grams have been written that do this analjrti- 
cally-^^-as 

Once superimposed, there are parts of the known 
structures that overlap very well, indicating that 
the structures are closely conserved in these regions 
(see bold areas in Fig. lb). We call these portions 
"structtirally conserved regions," or SCRs, and ex- 
pect that they will usually remain conserved in all 
members of this homologous family. These regions 
are usually composed of the secondary structure el- 
ements, the immediate active site, and other essen- 
tial structural framework residues of the molecule. 
Between these conserved elements are highly vari- 
able stretches that differ significantly from one 
member of the family to the next. These are called 
"variable regions," or VRs. They are almost always 
loops that lie on the external surface of the protein, 
and they contain all the additions and deletions be- 
tween different protein sequences. 

After the structures have been overlapped in 
three dimensions and parsed into SCRs and VRs, the 
next step is to align their amino acid sequences. For 
purposes of comparative modeling, the sequence 
alignment is done differently from the usual 
methods.^® Instead of relying on criteria such as 
amino acid identity or homology, we use strictly the 
three-dimensional overlap as the criterion. When 
the a-carbons of the respective residues in the over- 
lapped protein structures occupy the same place in 
three-dimensional space, then the residues are cor- 
responded in the sequence alignment (Fig. 2). Thus 
the alignment is primarily concerned with the 
SCRs, since these are the portions of the structures 
that are basically the same in all protein menibers of 
the family. The alignment in the VRs is often arbi- 
trary, unless two or more structures have similar 
conformations in a particular VR (for an example of 
this see the VR in the upper left of the structures 
fk' (sequence C-D-F-A) and «C" (sequence C-R-Y-V) 
m;Figs. 1 and 2). 

yj^The, resultant sequence alignment is then scruti- 




TABLE 11. Sequences of the Serine Proteases* 



COMPARATIVE MODELING TECHNIQUES 



Protein 



Chymotrypsin 

Trypsin 

Elastase 

Kallikrein 

Mast cell protease 

S. griseus trypsin 

Tonin 

Haptoglobin heavy chain 
Protein Z 
Protein C 

Nerve growth factor a chain 
Nerve growth factor 7 chain 
Blood clotting factor VII 
Blood clotting factor IX 
Blood clotting factor X 
Blood clotting factor XI 
Blood clotting factor XII 
Plasmin 

Apolipoprotein A 

Tissue plasminogen activator 

Urokinase 

Thrombin 

Complement factor B 
Complement factor 2 
Complement factor D 
Complement factor IR 
Complement factor IS 
Adipocyte protease 
Cathepsin G 
T-cell serine protease 
Hannuka factor 

Cytotoxic T lymphocyte protease 

Batroxobin 

EOF binding protein 

Drosophila snake locus 



Code Sourc e 

Bovine 
Bovine 
Porcine 
Porcine 
Rat 

S. griseiLs 
Rat 

Human 
Bovine 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Human 
Mouse 
Human 
Mouse 
Human 
Mouse 
Snake 
Mouse 
Drosophila 



CHTt 
TRPt 
ELAi 
KALt 
MCPt 
SGTt 
TONt 
HPH 
PRZ 
PRC 
NGA 
NGG 
VII 
FIX 
FAX 
FXI 
XII 
PLM 
ALP 
TPA 
UKH 
THR 
CFB 
CF2 
CFD 
CIR 
CIS 
ASP 
CAG 
TCL 
HAF 
CTL 
BTX 
EBP 
DSN 



*The sequences are taken from the sequence database » 
Tl-or those proteins for which a three^iimensional structure is 
tTe datwV'^'''"" "^^'^ """"^ Brookhaven struc- 



nized carefully to identify strong stretches or pat- 
terns of sequence homology that are characteristic 
of each SCR in the structure. For example, the se- 
quence "L-S/T-V-^-I-^," where tt is a charged or po- 
lar residue, is conserved in the second SCR in Figure 
2 Similarly, «G-LA» is found in all members at the 
third SCR. However, in the fourth SCR, the situa- 
tion IS less clear. A proline usually appears in the 
second position but may be replaced by a glycine 
Most SCRs can be identified by a characteristic (not 
necessarily contiguous) sequence pattern. In some 
cases, there is a pattern of hydrophobic or hydro- 
phihc residues rather than specific side chains. 
Careful analysis of the sequences and the corre- 
spondmg three-dimensional structures will usually 
allow some pattern to be discerned. 

The identification of these sequence homology 
patterns is a crucial step in the modeling. As is 
shown below, it is essential for the correct alignment 
of a new" sequence. If the characteristic sequence 
pattern is not present, then it is not clear how to 
align the "new" sequence and, consequently, how to 




construct the model structure. This is one 0^^^ 
r^n^v^' -nxparative modeling 
ren% depends so heavily on the retention bfSS 
quence homology. The above steps are perfbni^ 

only as new experimental structures and their se- 
quences become available. 

Construction of a "New" Structure 

tZ\T """"^ 'T^^ ^ ^^^^^ modeling a "new" pro- 
tein sequence of interest that is a clear member of 
this homologous family. The first step is to align this 
new sequence to the sequences of the SCRs using 
the previously identified characteristic sequence ho- 
mology patterns. It should be possible to align an 
appropnate portion of the "new" sequence to each 
f Thus, the "new" sequence 

L-i-V-K-I-E fits the pattern L-S/T-V-Tr-I-rr de- 
scribed above and can be aligned to the SCR at po- 
sitions 7-12 in Figure 2 and similarly for the G-I-A 
sequence of the SCR at 16-18. There will occasion- 
ally be more than one possible alignment for a par- 
ticu ar SCR. especially if the sequence pattern is a 
weak one. Such an example might occur in the SCR 
corresponding to residues 1-2 (Fig. 2). In such cases, 
the different possible alignments would have to be 
considered. As each portion of the "new*' sequence is 
allied to an SCR, the rest of the positions in that 
SCR are filled with the adjacent "new" sequence 
without permitting any additions or deletions 
within the SCR.t 

The remaining residues make up the VRs. The 
"new" sequence in each VR is examined to see if it 
corresponds in length and residue character to one of 
the VRs of the known structures. For example, the 
third VR In the "new" protein in Figures 1 and 2, 
residues 13-15, has the same sequence length and 
character as this VR in protein "B" (A-W-D-S-L in 
"new" vs. A-W-N-T-M in «B") and therefore has been 
aligned to it. On the . other hand, none of the known 
structures has a sequence that corresponds to resi- 
dues 3-6 of the "new" protein in residue length. 

Once the "new" sequence is aligned, as shown in 
Figure 2, the model building can begin. For the 
SCRs, usually the main chain coordinates fi-om any 
one of the known structures can be taken (Fig. Ic). 
The side chains are mutated to those of the "new" 
sequence wherever necessary. In our implementa- 
tion, the Xi angle for the "new" side chain is chosen 
to maximize overlap of this side chain on the old one. 
Further side chain x angles are not fitted at this 



1" An exception would occur if there are too few residues be- 
tween this SCR and the neighboring SCR. In that case, the end 
residues of an SCR would be left with a deletion. For examples, 
see position 96 in PLM and ALP in Figure 4. 
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Rgs. 1a-c. Leigend appears on page 321. 




Other 




re^<^ ^Lff'' ^^^^ *^ to the oext (VRs. 

ofSrSSifS^^^'^^**'^*^ c: Steps m me a>nstruct«>n 
if^Z^t^^"^ "^e* stfucnjf« are presented »n mts and 

f^i^c^LiT,^ ^> ^ constructed from me 

^v^^u ^ya^^ oisnyof^oithe known strucn^. 5irx» 
y are all very cios^ lo fhe wme T^e s*de cnams are mutated 



to &>e "f>ew" sequence where necessary, d: The various VR con- 
fofmatjons found in the known structures are considered for each 
VR of me "new" protein. {Compare the respective VRs shown 
here wrm those m a and b.) The ones mat do not fit are rejected, 
as shown by the crossed arrows. The most suitable is seiected in 
each case. In some cases, other conformational search"^-'*^ or 
energetics'*® methods must be employed, since no suitable con- 
format«>n can t>e found for mat VR among the known structures 
(see VR in upper tefl corner wim the sequence C-F-N-L-Q for an 
example in which a different, new conformation is necessary), e: 
The composfte structure showing me source of the respective 
"spare parts" selected for me model structure. 



# : 

A: 
8: 

C: 

NEW: 



jA 7}-: : - 



' J OiA K N T MjC 2 AlV - - A|E Gi 



V 



V r> I KjA WD 5 L.G I M 



js - - 



Fig. 2, Sequence alignment for me set of proteins in me sche- 
matic homologous family corresponding to Rgure 1, The align- 
ment is performed based solely on the overiap of me mree-dirnen- 
sionaJ structures and not by sequence alignment memods. the 
boxes delineate me SCRs as determined from the mree-dimeh^ 
sional stmcture overiap. The sequence of a "new" protein-JS 
aligned based on the characteristic patterns of sequence; hb^ol- 
ogy fourid for the known structures and meir sequen<^ (see text). 
The singie letter amino acid codes are given in Figure 4. . . :. 
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stage, because the nature of the different side chains 
usually diverges beyond this point. An automatic 
angle-scanning energy-optimizing routine is usually 
not employed, because we find that it often selects 
an unsuitable conformation based on minor or insig- 
nificant differences in energy. However, incorpora- 
tion of recently compiled and reported rotamer 
libraries'*" may improve our ability to automate the 
side chain conformation selection process. 

For those VRs for which an appropriate example 
structure can be selected from among the known 
structures during the alignment process described 
above, the main chain fragment is taken directly 
from the VR of that known structure, and again the 
side chains are mutated to fit the "new" sequence 
(see examples in Fig. Id and e). Remaining VRs, 
such as the VR at positions 3—6 in Figures 1 and 2, 
have no good example to build upon among the 
known structures. Therefore, they have to be con- 
structed by more complex conformation search 
methods'^^"'^^ and energetic considerations.*®'"**^ We 
have also used the protein database as a source of 
tentative starting conformations for loops of this 
type. The two ends of the particular loop that lie in 
the adjacent SCRs are taken, and the Brookhaven 
Protein Structure Database® is searched for struc- 
tural fragments that closely match the conformation 
of these ends (specifically, the a-carbon positions) 
and have the appropriate number of residues be- 
tween the ends. Thus, for example, in the VR at 
positions 3— 6 of the "new" sequence, the structures 
of the two ends corresponding to residues "A-D-A" 
(positions -1 to 2) and "L-T-V- . . ." (positions 7 to 
9. . .) are selected and typically ten structures that 
match these two ends best Qowest root mean square 
(rms) deviations] and that have five residues in be- 
tween are chosen and examined. Obviously, any 
such conformation that would collide with the re- 
mainder of the constructed "new" structure is elim- 
inated. Similarly, if the conformation does not pack 
well or buries a single charged group, it is rejected. 
In this way, one or a small number of initial confor- 
mations can be selected for this loop. Clearly this 
latter method does not produce an exhaustive list; 
however, it does provide loop conformations that are 
known to appear in other proteins. This method has 
previously been described by Kraulis and Jones*® 
for constructing protein structures from fragments 
using nuclear magnetic resonance (NMR) data or 
using crystallographically generated electron den- 
sity maps. 

It is important to emphasize that it is the initial 
spatial alignment of the three-dimensional struc- 
tures performed above (Fig. lb) that allows this con- 
venient clipping of fragments or "Spare parts" from 
the various known structures; to construct a compos- ; 
ite model structure of the "new*' protiein; Thus over- ^ 
^lapi^ing 6^^ structirres.: is essent^^^^^ 
steps: for the original.' aligninentkof ^the. sequences 



and for clipping together fragments in the subse- 
quent construction of each "new" model protein 
structure. Note that, the more experimentally 
known structures available, the more likely the 
boimdaries of the SCRs will be well-defined and the 
greater the likelihood that a suitable example 
known structure, i.e., spare part, will be found for 
each "new" VR among the known structures. 

Because the different portions of the "new" con- 
structed model structure arise from quite different 
sources (Fig. le), we can assign appropriate qualita- 
tive reliability confidence levels to the respective 
parts of the structure. Clearly, the conformations in 
the SCRs are the most reliable, especially when sev- 
eral known structures are available and have been 
compared in detail. This is true only if the respective 
SCRs of the "new" sequence retain the characteristic 
homologous sequence patterns in those SCRs. When 
they do not, and this is illustrated below, it should 
be regarded as a red flag warning that something 
different may be happening in this region of the par- 
ticular "new" protein. In such cases, great care must 
be taken in constructing the "new" structure and the 
confidence level in this portion of the molecule re- 
duced appropriately. For the VRs, when a good 
model structure appears for the respective VR 
among the known structures, then a fair confidence 
level can be assigned. This level is not as high as in 
the SCRs because of the greater inherent variability 
of these loop regions. The actual confidence level for 
each such VR would depend on the size of the loop 
(larger loops have more degrees of freedom and are 
therefore less reliably determined), how good the se- 
quence homology is to the known structure selected 
to model this VR, and how well the chosen confor- 
mation fits and packs onto the rest of the "new" 
structure. The confidence is lowest for those VRs 
with no good model loop among the known struc- 
tures. When conformational methods or database 
search techniques have to be employed, ex- 
perience^®"^® shows that we are not able to predict 
the conformation reliably in many of these cases. 

The resultant structure is then examined in detail 
to identify possible serious errors. For example, the 
inner cores of the structure are checked to be sure 
that no inappropriate charges have been buried. The 
modeling method makes it unlikely that there are 
serious steric contacts between main chain atoms in 
the model. This is because the main chain coordi- 
nates in the SCRs were taken from the known struc- 
tures. Only if the bad contact was present in the 
parent known structures will it be found in the 
model. In the VRs, absence of main chain overlap 
was a major criterion for acceptance of a conforma- 
tion. Therefore, in practice, bad steric contacts can 
usually be relieved by side chain rotations; occasion- 
■ally by having to select an alternative conformation 
for a VRV The moder structure can then be intro- 
fduced into' an energy-minimization program'*®*'*'^ to 




Fig. 3. a-Carbon plots showing the superposition of the seven 
experimentally known serine protease structures (see Table I), 
Note the SCRs, where the seven stnjctures are virtually the same, 
and the VRs, where the structural divergence is considerable. 



The key to the plots is as follows: CHT, green; TRP, solid red 
ELA, solid blue; KAL, solid cyan; MCP, dotted green; TON. dottec 
red;SGT, dotted blue. 



relieve any remaining minor steric contacts and to 
optimize the bond angles and torsional angles. Typ- 
ically, we begin with restrained or template-forced"^^ 
minimization, forcing the atomic coordinates to re- 
main close to their initial positions. This allows bad 
contacts to relax without introducing large artifac- 
tual distortions into the structure. After several 
hundred cycles, the constraints can be progressively 
removed and the energy of the molecule minimized. 

Experience with comparative molecidar mod- 
eling^^-^3 has shown that these models are never 
completely accurate, especially in the conformations 
selected for the VRs. Therefore, any uses of the 
model structure and predictions that are made based 
on it should take into consideration this fact and the 
local reliability confidence levels discussed above. 
We always regard such a model as a working hy- 
pothesis; as new data emerge, whether from struc- 
tural, spectroscopic, or biochemical sources, the 
model is modified, corrected, improved, and refined 
to fit the new data. 

RESULTS AND DISCUSSION 
Initial Alignment of Known Serine 
Protease Structures 

Our analysis of the serine proteases began some 
years ago with the superposition of the then three 
known structures of chymotrypsin, trypsin, and 
elastase.^-* Since then, coordinates for four more 



structures have appeared in the literature and data 
bank®: kallikrein,^'^ mast cell protease,^ Streptomy- 
ces griseus trypsin-like protease,^^ and tonin.^^ The 
overlap of these additional structures was performed 
by first superimposing the added experimental 
structure onto the rest using the a-carbon positions 
of residues 195, 57, and 102.t The superposition was 
refined, when necessary, using repeated cycles of 
least squares fit on related residue a-carbon posi- 
tions between the protein to be fit and chymotryp- 
sin. Only residues whose positions were conserved in 
the two structiures were fitted in the calculation; 
that is, any residue was eliminated from the fitting 
if the deviation in a-carbon positions was more than 
2 A on the first cycle and then more than 1.5 A on 
subsequent cycles. Figure 3 shows the final overlap 
obtained for these seven structures. Virtually all the 
SCR positions identified in the previous analysis of 
the first three proteins* remain structurally con- 
served. The repertoire of conformations, i.e. "spare 
parts," found for each of the VRs is increased with 
the larger number of structures, as expected (Table 

ni). 



tThe residue numbering given for the serine proteases fol- 
lows that of the chymotrypsinogen molecule throughout this 
paper, in the text, in the tables, and in the. figures. 
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TABLE III. Sequence Lengths of the Different VRs in the Serine Proteases 

VR 



VR 


35-41 


59-62 


72-80 


97-101 


125-133 


146-151 


166-179 


185-189 


203-206 


217-224 


Orll 


7 


4 


9 


5 


9 


4 


14 


3 


4 


8 


TKr* 


5 


3 


9 


5 


7 


6 


14 


5 


0 


8 


ELA 


10 


5 


9 


7 


9 


5 


16 


4 


4 


•10 


XT' AT 


6 


3 


9 


9 


7 


8 


14 


5 


0 


9 


ML-r 


10 


3 


9 


5 


9 


5 


13 


5 


0 


6 


SGT 


2 


8 


7 


3 


7 


5 


15 


6 


5 


8 


TON 


4 


3 


9 


16 


7 


6 


14 


5 


0 


9 


HPH 


6 


16 


0 


1 


8 


4 


31 


5 


6 


7 


PRZ 


7 


1 


4 


5 


11 


2 


13 


3 


4 


6 


PRC 


7 


5 


9 


5 


13 


10 


14 


5 


4 


8 


NGA 


6 


3 


9 


16 


7 


6 


14 


5 


0 


9 


NGG 


6 


3 


9 


12 


7 


6 


14 


5 


0 


9 


VII 


Q 


g 


9 


5 


1 9 


5 


1 Q 


c 

Q 


A 
*± 


g 


FIX 


6 


5 


9 


7 


11 


5 


14 


5 


4 


8 


FAX 


7 


5 


8 


5 


11 


5 


14 


5 


4 


8 


FXI 


9 


8 


9 


6 


9 


5 


15 


5 


3 


8 


XII 
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Based on the above superposition of structures, 
the seven sequences were aligned as described above 
(Fig. 4). Note the wide variation in the length of the 
different sequences in each of the VRs (Table III). 
Using this sequence alignment, we attempted to 
identify the critical conserved stretches of amino 
acid sequence. These are characteristic of each SCR 
and are essential for aligning "new" sequences to 
the known ones, which is the first step in modeling. 
. The patterns that were chosen are noted in Figure 4 
at' the bottom of the sequence list in the row labeled 
I'jCpN.'* In some cases, there was almost complete 
. conservation of a single amino acid. In other cases, 
a hydrophobic side chain (typical of an internal 
^Mfe5'^.;^.S -position), seems to be required. 
;,'':^,H'!;*: "\V?';Becaiise the identification of these stretches is 
j.'') /Siich an important part of the modeling process, it is 
i^^^4|v.w:6rthwhile describing in more detail how this: is 
J^^^SSj^rfbraied: Iri most cases, it. is straightforward;^ the 
ilQS^i^&iasei^ed' seqiiehce is readily apparent.^ i^or,. exam-; 
■^"^^ll^l^^tjh^AcharaV^^ UV-I/V-G-G at positionscie^^. . 

^tev:. ■■■ ■ ... 

^^i;. . ... 



19 (see Fig. 4) defines the first SCR very clearly. 
Similarly, the remarkably conserved S/T/A-G-W-G 
sequence of residues 139-142 is highly characteris- 
tic of this SCR. On the other hand, the SCRs at 63- 
71 and 81-96 have no such clearly conserved se- 
quence, and it is often very difficult to align a "neV* 
sequence to the known ones in these SCRs unambig- 
uously. 

To understand better why this latter group of 
SCRs has so little sequence homology, the structures 
of these two regions were examined in detail (Fig. 5). 
It is immediately apparent that these SCRs lie on 
the surface of the molecule, and virtually all the 
residues point out into solvent and thus can and do 
vary with little consequent effect on the rest of the 
structure. However, in each SCR, there are some 
side chains that point into the. interior of the mole- 
cule? and thus are more , conserved; In the SCR at 
63:^71; .1 this includes ■ the: ,Val |;(or';. comparable ali- 
. phatic residue) at 66, and . the aliphatic side chain at 
68,TTjie!,Gly at position 69.iappears to be almost com- 
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Rg. 5. The location is shown for some of the SCRs that are 
difftcuft to align in the serine proteases because of only weak 
characteristic patterns of hotnotogy. This is usually due to the 
SCR lying on the surface of the molecule. This figure presents the 
Ca plot of CI-fT (dotted lines), with the location shown of the 
respective SCRs that are difficult to align: 63-71 and 81-96 (solid 
lines). Note that these SCRs lie on the surface of the molecule and 



that most of their side chains point out into solvent and thus are 
free to vary more rapidly than internal residues. In the SCR at 
63-71, conserved residues 66 and 68 are the only internal ones 
aruj are always hydrophobic; 69 is always a Gly, because it forms 
a turn. In the SCRs at 81—96, internal hydrophobic residues are 
83, 85, 89, and 94. Partial conservation is observed at positions 
91 and 92, with a His-Pro. 



pletely conserved, probably because it forms a turn 
and adopts a conformation (<}> = —60°, »j> — —25**) 
that is permitted only for Gly (Fig. 5). Similarly, 
small hydrophobic side chaiins occur at positions 83 
and 85 in the SCR at residues 81-86 (Fig, 5). More 
varied but always hydrophobic residues appear at 
89, and an aromatic one is typical at 94, both of 
which point into the molecule, in addition, there is 
also a tendency for a His-Pro at positions 91—92 (Fig. 
5)i however, a number of sequences have only one or 
none of these two residues (Fig. 4). 



CNO 55 60 65 70 75 80 CNO 

> < ><-> <- 

CHT AAHCGVT TSDWVAGEFDQGSSSE-KIQK CHT 

TRP AAHCUKS — GIQVRLGQDNIMV-VEGNQQF TRP 

ELA AAHCVDRE LTFRWVGEHNLNQ-NNGTEQY ELA 

KAL AAHCK NDHYEVWLGRHNLFE-NENTAQF KAL 

MCP AAHCKG REITVILGAHDVRK-AESTQQK MCP 

SGT AAHCVSGS GNNTSITATGGWDLQS GAAVK SGT 

TQM AAHCYSH HYOVT.LGPt>IHT.FK-DEPFAOR IQH 

HPHoid TAKMLFLN"=— HSENATAKDIA-PTIiTLY—VGKKQL HPHold 

HPHnew TAKHLFLNHSENATAKDIAPTLTLYVCKK"-""-— QL HPHnew 

> < ><-> <- 

CON AAHC V XG CON 

Rg. 6. The old^-* and new sequence alignments for HPH, rel- 
ative to throse of the known three-dimensional stnjciures, for the 
region around the SCR 63-71 . Note that six residues have been 
moved from the VR after position 71 to the VR prior to position 63. 
Nomendature and latieling are as in Figure 4. 



: ; Rg. 4. Sequence alignment for many of the known serine pro- 
V - tease sequences. The source of these various sequences and the 
. ^definitions of the protein name codes used In this Figure are given 
. . - in Table II. The first seven proteins (above the line) are those with 
; ^ known experimental three-dimensional structures (see Table I). 

> Vj . VTheir sequences have been aligned based on the superposition of 
r : Tjthe three-dimensional stmctures as described in the text. The 

• \=;, kemalningj proteins were aligned using the characteristic se- 

i-'quence matching patterns (see text). When the sequence align- 
; ^:^ment ^s uricert^ the respecUve sequence Is shown in italics. The 
•^^,V^'i^^UPACrlUB convention standard single letter amirra acid code. 
') --'Tused In-thls figure^ is as follows: A = Ala, C = Cys; D - Asp, E 
Glu,.T:^ Phe.G = Gly.H = His, I = lie. K = Lys, L = Leu, 

• • i = :M^t;:N = Asn, P = Pro. Q = Gin, R ,= Arg, S = Ser, T = . 
^\H^^^I#V^ V^= 'Val. W = Trp, Y = Tyr. An asterisk indicates that the. 

;<ii|v^t^s«?Oefee>-OT^ past the. last residue: shown in the figure^ 
C.^^^f^PositiG^ns of relative dele^^ the sequences aro denoted by a 

■ ■■■ - ■ ■ ■ 



dash in the known structures and by a double dash in the se- 
quences aligned using the characteristic homology patterns. The 
line labeled "CNO" gives the chymotrypsinogen residue nunrUjer- 
Ing used throughout this paper. The symbols delineate the 
SCRs. The "CON" line lists the conserved, characteristic se- 
quence patterns used to align the "new" sequences (see text). 
They are coded as follows: upper case; almost completely con- 
served side, chain; lower case, high frequency of this amino acid 
at this.' position; X, nonpolairreslduesi typically A, V, L, 1, M; ir, 
denotes -a^ polar residue, typically ^S.^T,'^ Qv N plus the charged 
resrdues;!p,rdisnotes.S or. T,v'(Fnay:te^^^^ by A occasion- 

ally)'; <t>'.sderiotesrarbmatic, usually^iY- dr. Fi sometimes W; +, de- 
notes a positive residue such; as-ft; Ki or H; r-, denotes a negative 
residue such as D or E; The line labeled. "S-S" gives the position 
:'.\ otthe^ha(i;cystiriei,readue:i1hM fe,)^^^ to the half 

>;:^Uivc^tne,a^^ positions 
^v;^' where^itds nold bridge is formed. 
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Rg. 7. Ca plots for the old (dotted Ifnes) and new (solid lines) 
fl^nr^Trl ^^.^ ^^'^ °" respective alignments shown fn 
iTy^ o J '^'l® '^"S^f ^PP®^'' around the SCR at positions 
63-71. Because the VR from 59 to 62 in the new HPH alignment 
IS SIX residues longer than in the previous alignment (Fig 6) and 
at least eight residues longer than any of the currently known 
s^njctures .n this loop (see Figs. 4 and 6, Table III), the confor- 
mation that has been drawn for this VR (dashed lines) is arbitrary 
Lf^^Vif-^T included only to give a sense of the relatively large 
size of this loop. The difficulties associated with finding the correct 



Th^ wS Z '^''9® additions are described fn the text 
Il^ortPr in t^fn''''^^^^".^° '^.^'^^'^ ^^f^^ated. seven residues 
t\^?Il!!' Compare the conformations of the 

known staic ures for the VRs at 59-62 and 72-80 (Fig 3) wrth 
mese loops in HPH Whereas the Ca positions (and n^l n chS n 
coordinates) for residues 63-71 are the same In this figure for ?>e 
old and new structures of HPH. the side chains at each residue 
position are. of course, different in the two alignments (Rq 6) 
giving different overall stmctures for this part of the molecule as 



Alignment of "New" Serine 
Protease Sequences 

Identification of the characteristic sequence con- 
servation patterns for each SCR allows the process 
of aligning "new" sequences to begin. Figure 4 
shows the alignment of a large number of serine 
protease sequences to those of the known structures. 
The proteins included are taken from a very wide 
variety of functions and species. It is worth noting 
that three of these proteins, although clearly homol- 
ogous to the rest and therefore full members of the 
femily, are no longer functional serine proteases. 
These are the heavy chain of haptoglobin (HPH) 
protein Z (PRZ), and the a-chain of nerve growth 
factor CNGA). For all these proteins, the character- 
^ic pattern for each SCR was located and matched. 
Then the remaining positions in that SCR were 
tilled, without permitting any additions or deletions 
One of the results that emerges immediately from 
this alignment (Fig. 4) is the strong reinforcement of 
the characteristic patterns in all these "new" se- 
quences. The same basic patterns appear in virtu- 
ally every protein, yet, in almost every pattern, 
there are individual protein sequences that have 
exceptions. Usually the exceptions are minor devia- 
tions from the theme of the conservation. For exam- 
ple the chosen four-residue stretch from 139 to 142, 
b/r/A-G-W-G, is present in 26 of 35 protein se- 
quences reported in Figure 4. The Trp is replaced by 
a Phe or a.Tyr, both of which are typically accept- 



able replacements for a Trp, in four other proteins. 
In the same way, the characteristic sequence C-G-G 
almost always appears at positions 42-44, but some- 
times one of the Glys becomes an Ala, and occasion- 
ally the Cys is replaced by an alternative small side 
chain. It is clear that, with perhaps the occasional 
exception of the difficult SCRs between 63-71 and 
81—96 discussed above, all the sequences shown can 
be aligned relatively trivially and unambiguously to 
the respective characteristic conserved sequences in 
the SCRs (Fig. 4). 

The close examination of the SCRs at 63-71 and 
81-96 described in the previous section (Fig. 5), to- 
gether with the strong reinforcement of the charac- 
teristic sequence patterns in these SCRs (Fig. 4), has 
led to the realignment of the HPH sequence in this 
region. This results in a large change in the three- 
dimensional model structure for this protein from 
that previously proposed.^*^ The old alignment (Fig. 
6) was forced by the need to avoid placing the 
charged Asp 65 in position 66, where it would be 
buried in the hydrophobic core of the N-terminal 
3-barrel. Pro was placed at position 69 as a replace- 
ment for a Gly that would permit a turn. Realization 
that the aliphatic residues at positions 66 and 68 
and the almost complete conservation of a Gly. at 
position 69 (Figs. 4 and 5) make up the characterisr. 
tic sequence pattern for this SCR led to the realigni 
ment of the HPH sequence relative to those :}bf\fiti# 
known structures as shown in Figure 6.. As a result; 
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Fig. 8. Ca plot for CHT is shown indicating the sites of VRs 
with large additions (dashed tines) and with large deletions (dotted 
lines) in CFB. The unusual residues described in the text are 
represented by solid lines. These residues include 16-20, 43, 
139-142, and 154-156. As can be seen, the unusual residues 
are localized in one region of the molecule around the normal site 



of the N-terminus of the protease peptide chain, residue 16. Res- 
idues 1 6 and 1 7 should not be in their usual positions, because 1 6 
is not N-terminal in CFB or in CF2. Consequently, they really 
should be placed In the approximate position expected from our 
knowledge of the zymogen structures'"^ where the deavage at 
16 has not yet occurred. 



of the new alignment, the HPH sequence has a very 
large addition in the VR at 59—62 ajid a very short 
loop in the VR at 72—80. Now, however, the require- 
ments of aliphatic residues at 66 and 68 and Gly at 
position 69 are completely satisfied. The difference 
between the old and new structures for HPH is il- 
lustrated in Figure 7, which demonstrates clearly 
that changes in the alignment can result in large 
changes to the derived structures. 

From an evolutionary perspective, deviation from 
the typical pattern of sequence conservation occiirs 
most often in those members of the family that are 
no longer serine proteases in fvmction. For example, 
HPH and PRZ are no longer serine proteases in 
function. Both have lost the active site Ser 195 (it is 
: an Ala), and His 57 is also replaced by a Lys in HPH. 
H.. Nevertheless, both proteins retain most of the char- 
■ V :■ acteristic sequence homology patterns, yet they de- 
• V: / vyia'te;,inore than the typical sequence. Examples of 
. ;V this;;a^ immediate area around the Ser 195 or 
. -li r tihterlsfek ^it 63:^71. The Streptomyces griseus trypsin- 
V /i K^^^^^^^ a weaker homology pattern 

. : ' ■ ^aLn;|fief testv^^ as has been previously noted.'* Even 
^ ;i' -;^;;>lu^^ are: the other bacterial serine 

; " X J;,'/' ]pfdt^jasfes.?^,.These.latter are not considered iii this 

y.^>;\t^a^^^ considerably great6r-dijEiicuJty : 

/:i>5i^$vs;iun^^ correctly I by;- - thiese'. ■ ) 

^'i^S$^^mMik^^ oiirpage.-318).= ^:..yk-^ ■ 



Analysis of the Variable Regions 

Once the sequences that correspond to the SCRs 
are identified and aligned, attention can be focussed 
on the remaining residues forming the VRs. Based 
on the alignment in Figure 4, the sequence lengths 
deduced for the various VRs are summarized in Ta- 
ble ni. It is clear from Table III that an amazing 
degree of variation occurs in the length of these 
loops and therefore in their conformations (Fig. 3). 

In previous studies,* we classified the modeling of 
the VRs into five diflferent classes. 1) A known struc- 
ture has the same-length VR and the same residue 
character. In this case, a model structure is con- 
structed for the "new" VR by clipping in the respec- 
tive VR from the known structure. In such cases, the 
confidence level in this part of the structure can be 
reasonably good, especially if it is a short VR. One 
example occurs among the known structures of the 
serine proteases at residues 97-101, where the VRs 
of chjonotrypsin (CHT), trypsin (TRP), and the mast 
cell protease (MCP> arc all : the same length and 
ha.ve the same ^conformation? (see this loop in Fig. 3). 
2) Different lengths ; app6ar^in{^ among the 

known stifucturesy,bu^^ 

tural^ mo tif i ; e; gyi' BL' . pTbend^ai?necting^:fewo i antiparal r 
lei ;^TSt^ands^?^^ 

vtui^ls^^^ be 
^villustraf^bjr^^^ainining^^ positions 35—41 
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the DSN protein in the SCR from 63 to 71. The usual 
Gly residue that invariably appears at position 69 is 
absent. This leads to an ambiguity regarding how to 
align the sequence. However, at least one of the pos- 
sibilities, the one shown in Figure 4, is likely to fit 
into the usual serine protease structure at this po- 
sition with only minor deviations. 

On the other hand, there are examples of larger 
changes, including the complete absence of one or 
more characteristic sequence patterns. This can be 
illustrated by the SCR at 134-145. This SCR has the 
highly conserved, very characteristic S/T/A-G-W-G 
sequence at 139-142, as noted above. Virtually ev- 
ery serine protease has this sequence, with a rare 
change from Trp to Tyr or Phe and a very rare 
change of the second Gly to something else (see PRZ 
for an example of divergence as discussed in the pre- 
vious paragraphs). Note, however, the sequence for 
this SCR in CFB and CF2 (Fig. 4). In these cases, the 
characteristic residues are completely missing. 
Clearly something is happening to the structure 
here that is unusual. The absence of these residues 
creates a serious problem in properly aligning these 
two sequences in this SCR. There is no obvioxis pat- 
tern of residues that would correspond to the S/T/ 
A-G-W-G sequence at 139-142 among the residues 
that are in this region of the sequence. It is possible, 
in principle, to place the sequence A-L/H-F-V at 
these positions, thereby replacing the conserved Trp 
with a Phe. However, this leaves internal positions 
136 and especially 138 with charged residues such 
as Asp and Lys, where nonpolar residues are char- 
acteristic and are called for by the structure. The 
sequences have been tentatively aligned as shown in 
Figure 4 to satisfy the requirements that internal 
residues be nonpolar, or at least uncharged. The un- 
certain alignments are shown in italics in Figure 4. 

If we look further at these two sequences, CFB 
and CF2, we find another anomalous site. The clas- 
sic, characteristic pattern for positions 42-44 is C- 
G-G. Examination of the sequences in Figure 4 
shows that substitution of one of the Gly residues by 
Ala occurs occasionally. The startling introduction 
of a much larger Met side chain at position 43 in 
CFB and an even larger and charged Arg in CF2 
indicates that once again something strange is hap- 
pening in these proteins. Examination of the loca- 
tion of -these changes in the serine protease three- 
-dimensional structure (Fig. 8) shows that they are 
immediately adjacent to each other in the structure 
even though they are qxiite distant in the sequence. 
:- :j;?VA further anomalous sequence occurs at positions 
16-19 of these two sequences^ Activation of a serine 
piTotease from, its zymogen i protein in vo Ives the 
. cleavage of the peptide chain, prior to residue 16 to 
■ generate a free; amino terminus at this point ia the 
;,m'olecule,^^?® The earliest , serine protease ^crystal 
^^smictui'e solutions^'* showed that this amino: termiT : 
ifeu^j foiined; an!;; internal: salk bridge : with, the side ' 



chain of Asp 194. The homology observed in the res- 
idues at positions 16-19 is likely to be diagnostic of 
this activation step in all these enzymes. Thus it is 
interesting that several of the sequences in Figure 4 
show major deviations from the characteristic se- 
quence pattern for these four residues. All the serine 
proteases have He or Val (or Leu) at positions 16 and 
17 and usually Gly at 18 and 19. Exceptions to this 
pattern are found in PRZ and NGA, neither of which 
is a true protease and thus presumably no longer 
retains an activation step. The two proteins CFB 
and CF2 are also striking in the absence of any ho- 
mology to the characteristic pattern for this stretch 
of residues or even between the two themselves. It is 
known that, unlike the other serine proteases, these 
two proteins are clipped some 200 residues N-ter- 
minal to residue 16.^^*^* Thus the activation mech- 
anism for these proteins must be different from that 
of the other serine proteases. Figure 8 shows that 
the normal location of the N- terminal residues 16 
and 17 is also close to the two unusual sequences 
described above for the three-dimensional structure. 
Thus significant violations of three spatially ad- 
jacent characteristic sequence patterns appear in 
these two proteins, coupled with ftmctional signifi- 
cance expressed in the absence of the normal mech- 
anism of zymogen activation. It is interesting to note 
that CFB and CF2 occupy parallel functional roles 
in the alternative and classical pathways of the com- 
plement system, respectively.®^ They form the pro- 
teolytic subimits of the respective C3 convertases in 
the two pathways and, upon binding of C3b sub- 
units, change their specificity to become C5 conver- 
tases. Taken together, these observations suggest 
that some unusual three-dimensional structure is 
occurring in this particular region of the molecule 
involving the activation site and the immediately 
adjacent active and specificity sites. What is difficult 
to predict or determine is just how different these 
structures actually are from the typical serine pro- 
tease structural theme in this region. The excellent 
preservation of the characteristic sequence homol- 
ogy patterns everywhere else in the CFB and CF2 
sequences (Fig. 4) suggests strongly that the re- 
maining parts of the molectde are close to their nor- 
mal conformation. However, there is no such reas- 
surance for this part of the molecule. 

As one begins to construct this part of the CFB 
molecule onto the typical SCR framework of the 
serine proteases, several problems are encountered 
very rapidly. The introduction of side chains at po- 
sitions 140 and 43: in place of Gly, in such close prox- 
imity, causes a steric collision of these two groups. 
In cohstrutting the N-terminal portion of the chain, 
tl^e conformation for residues 16-19 must be: taken 
£rpm: :the, zymogen - structure of chymotrypsin- 
• • residueis cannot reside in the 

J;V^^ji^^al N-tehniiial he 16 biridifig pbcket^'This leaves 
to "do with Asp a94 the other 



COMPARATIVE MODELING TECHNIQUES 



half of the normal salt bridge to the amino terminus 
at He 16. Both possible conformations of the aspar- 
tate, that of the zymogen and of the mature enzyme, 
must be considered based on zymogen crystal 
studies.^'^-''** Unfortunately, neither conformation 
seems possible without a significant change in the 
local CFB structure. In the zymogen conformation, 
the side chains of Met 43 and Ala 140 (both usually 
glycines) prevent the Asp 194 side chain from occu- 
pying its usual place. In the mature enzyme form, it 
is the side chain of Phe 142 (also typically a glycine) 
that is too close to permit the usual conformer (see 
Fig. 8). 

The resulting changes that need to be made to 
accommodate these unusual side chains are difficult 
to predict. Does the structure undergo a large num- 
ber of small changes that results in the accommoda- 
tion of the above groups in close to the normal con- 
formation, or is there a significant change in this 
portion of the molecule that gives a very different 
conformation around Asp 194 and residues 139- 
142? It is not possible to answer this question using 
the energetics methods currently available. Such 
problems may have to wait either for experimental 
structure determinations or for improved under- 
standing of protein structure and energetics. There- 
fore, it is particularly important that comparative 
modeling methods, as described in this paper, allow 
us to recognize such regions of the molecule that are 
likely to be particularly difficult to construct reli- 
ably. 
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